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Abstract 



In order to test the canonical quantization programme for general relativity 
'nT , we introduce a reduced model for a real sector of complexified Ashtekar gravity 

which captures important properties of the full theory. While it does not corre- 

^ ■ spond to a subset of Einstein's gravity it has the advantage that the programme 

^Q . of canonical quantization can be carried out completely and explicitly, both, via 

(--^ \ the reduced phase space approach or along the lines of the algebraic quantization 

^^ ' programme. 

—.^ ■ This model stands in close correspondence to the frequently treated cylindrically 

Q>^ I symmetric waves. 



In contrast to other models that have been looked at up to now in terms of the 
O^' new variables the reduced phase space is infinite dimensional while the scalar 

^ , constraint is genuinely bilinear in the momenta. 

The infinite number of Dirac observables can be expressed in compact and ex- 
plicit form in terms of the original phase space variables. 
^ , They turn out, as expected, to be non-local and form naturally a set of countable 

H ' cardinality. 

1 Introduction 

Ashtekar's variables ([|I|) simplify the algebraic structure of the constraint equations of 
general relativity so tremendously that one can solve various problems of classical and 
quantum gravity that were simply infeasible before in terms of the old ADM variables. 
In particular, a number of model systems could be solved to more extent or even com- 
pletely in these new variables. 

The model systems which could be solved completely up to now in terms of the new 
variables (e.g. [@]), to our knowledge, lack from not having in common at least one of 
the following features with general relativity : 
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1) the reduced phase space is infinite dimensional (i.e. the number of Dirac observ- 
ables) and 

2) the scalar constraint is bilinear in the momenta. 

Here we introduce a new model which is an algebraic plus Killing reduction of com- 
plexified Ashtekar gravity. In order to account for the fact that Ashtekar's connection 
obeys a non-trivial reality condition, we also impose reality conditions that bring us 
'as close as possible' to the cylindrically symmetric waves (0) which form a genuine 
subset of Einstein gravity whereas our model defines then a real sector of complexified 
gravity which is manifestly no subset of Einstein gravity. 

However, our model has both features 1) and 2) and we are thus able to attack one of 
the major problems of general relativity, namely to find a complete set of (preferredly 
canonically conjugate) Dirac observables for a theory with an infinite number of degrees 
of freedom and non-trivial dynamics. This is important in order to get some intuitive 
insights for the problem of non-perturbatively quantizing full general relativity. 
The results of this paper are as follows : 

In section 2 we define the present model and show that it is not quite possible 
to formulate the cylindrically symmetric waves (which means to project our model 
of complex gravity to the real sector corresponding to Einstein gravity) in terms of 
Ashtekar's variables after we carried out the reduction. We then impose reality condi- 
tions that bring us at least 'as close as possible' to the cylindrically symmetric case. 
In particular, the first- and the diagonal part of the second fundamental form are real 
while the off-diagonal part of the latter is purely imaginary. 

In section 3 we carry out the symplectic reduction (0]) of our model. 
That is, we are using the degenerate Hamilton- Jacobi method to solve the constraints. 
Quite surprisingly, it is straightforward to solve the typically very complicated scalar 
constraint before solving the diffeomorphism constraint. As usual for Hamiltonians or 
scalar constraints bilinear in the momenta, there are several 'sectors' of the constraint 
surface. One of them corresponds to degenerate metrics, the other to non-degenerate 
metrics (the latter defines the 'physical' sector of our model). 

The most complicated step is then to find a complete but not over complete set of 
diffeomorphism invariant observables. This turns out to be feasible and the resulting 
functionals of the basic phase space variables can be given in explicit and compact 
form (the observables given in ^ are rather complicated and intractable objects due 
to the fact that the mode system of Bessel functions is involved). 
As expected for a diffeomorphism invariant field theory, the Dirac observables come 
out to be naturally smeared and nonlocal. In contrast to non- diffeomorphism invariant 
field theories one does not have to employ a complete system of mode functions on the 
hypersurface to write the Dirac observables in this form. The reduced phase space thus 
has naturally a countable cardinality. We argue that this is always true for a spatially 
diffeomorphism invariant field theory. 

Section 4 is devoted to quantum theory. 
From the reduced phase space point of view this is already trivial since the phase space 
derived in section 3 can be cast into a form so that it has the structure of a cotangent 
space over a real, infinite dimensional vector space. 

The application of the algebraic quantization programme forces us to make use of func- 
tional integral techniques. 



The two resulting Hilbert spaces turn out to be unitarily equivalent. 
As to be expected for a field theory with constraints quadratic in the momenta, the 
final set of Dirac observables mix configuration and momentum variables with respect 
to the original polarization (here the Ashtekar polarization). 

While it is important to derive the Dirac observables, they are typically hard to inter- 
pret. Here the framework of deparametrization (0) helps to define 'evolving' Dirac 
observables which allow for a direct interpretation in terms of metric variables and to 
leave the limits set by the 'frozen' formalism. 

We extend this formalism to field theories and apply it to our model. We thus arrive 
at an interpretation of the theory although the meaning of the constants of motion 
(Dirac observables) remains veiled. 

Finally, in section 5, we discuss the relation to the loop representation (J^) of our 
model. 

First we define a set of loop variables that stand in bijection with the phase space 
variables underlying the connection representation. This implies that we can explicitly 
express constraints, symplectic form and Dirac observables in terms of loop variables. 
These expressions are quite complicated compared to those obtained for the connection 
representation and that allows for two possible conclusions : 

Either it seems to be doubtful whether the loop variables are tailored for the problem 
to construct the reduced phase space or 

Killing reduced models are unsuited to simulate the situation in the full theory as far 
as the loop representation is concerned. 

We summarize what has been learnt in dealing with the present model, valid for 
the full theory of Einstein gravity. 

2 Definition of the model 

We assume the reader to be familiar with the Ashtekar formulation of canonical gravity 
(|I[]). In the sequel a,b,c,.. will denote tensor indices and i,j,k,.. SO (3) indices and we 
apply the abstract index formalism. 

We impose the following three conditions on the phase space variables of the complex- 
ified Ashtekar formulation of gravity : 

1) Killing-reduction : all fields are Lie- annihilated by a set of two commuting, hyper- 
surface orthogonal Killing fields. 

Since the Killing fields commute they are surface forming. Since they are hypersurface 
orthogonal these surfaces lie orthogonal to a foliating vector field. Hence we find three 
globally independent coordinates x,y and z such that the Killing fields are given by 
dx, dy. Then all fields depend on z (and the time variable t) only. 

2) algebraic reduction : We set to zero the following components of the basic fields of 
the Ashtekar formulation 

Al = Al = Al = Al = El = El = El = El := . (2.1) 

Up to now we are on the phase space of the so-called Gowdy models ([§]). We will 
assume, for simplicity, the spatial topology of the initial data hypersurface to be that 
of the three-dimensional torus T^ in order to avoid technical difficulties associated with 
boundary conditions. 



One more condition brings us close to the axisymmetric spacetimes. 
3) algebraic reduction : The Killing fields are orthogonal. 
This last condition forces the densitized triad to be of the form 

E^ cos(/3) -Ev sin(/3) 
{E1)=\ E^sin(/5) E^cos(/5) | (2.2) 

Q E^ 

and in order to have as many configuration variables as momentum variables we require 
in analogy to (2.2) the Ashtekar connection to reduce to 

Ax cos(a) —Ay sin(a) 
{A\)=\ A^sin(a) AyCos{a) | (2.3) 

A, 

where all variables are generally complex up to now. 

As far as the triads are concerned, provided we restrict them to be purely real, we 
easily see that we have already a real diagonal 3-metric (compare the appendix). Let 
us see whether we can also implement the reality condition 

Ke{A\ - r^) = 0, r^ = spin-connection (2.4) 

which leads to the Einstein phase space. In order to do that we first have to compute 
the spin-connection. The calculation is deferred to an appendix. The result is given by 

/ -r,sin(/3) VyCos{l3) 
{V\)= r.cos(/3) r,sin(/5) I (2.5) 

V or, 

where 

r, = ^(±^j^)' r = — ^—{t^-E^y r, = -i3' . (2.6) 

2Ey ^ E^ 2E^ ^ Ey 

The problem is already obvious : both, F^ and A\ have five non vanishing components, 
giving rise to five reality conditions induced by (2.4), but these components are coor- 
dinatized respectively by only four independent coordinates (namely A^., Ay, Az, a and 
E^,Ey,E^,j3). The explicit form of these reality conditions is given by the require- 
ment that the following quantities (i-ft'* := iKabE'^/Jdet{q) = A^^ — T^^ where Kab is 
the extrinsic curvature) 

iK^ := Ax cos(a) + F^ sin(/9), iK'^ := A^ sin(Q;) — F^ cos(/9), 

iKy := —Ay sin(a) — Ty cos(/9), iKy := Ay cos(q;) — Ty sin(/3), 

iKz := Az - Tz (2.7) 

are imaginary. Equivalently 

tKi^Ef = Ey{Ax sin(a - /?) - FJ, tK^E^ = E^{-Ay sin(a - /?) - F,), 

iKlE"^ = E'^Ax cos(a - (3), iK^Ef = E^Ay cos(« - (3), 

iKiEt = E\Az-Vz) (2.8) 



are imaginary. 

We will now show that these are five functionally independent equations whence our 

reduction is incompatible with (a subset of) the real form of complexified gravity that 

corresponds to the Einstein phase space : 

We will immediately show that the Gauss constraint of this model is given by 

g = {E')' - sin(a - f3){A^E'' + AyE^) . 

Let 5 := a — [3. Then, using the Gauss constraint we derive that 

^''''^^^ ^ E-[A,cos{5)] + Ey[AyCos{5)] ^^'^^ 

which is weakly imaginary since the square brackets in the denominator are according 
to (2.8). Assuming that the imaginary part of 5 is finite, this implies unambiguously 
that 5 itself is imaginary : S = iC, where ^ is real. Then cos{6) = ch(^) is real whence 
again from (2.8) we derive that A^, Ay are imaginary. Finally we have j := A^ + a' = 
{Az — r^) + 5' = {Az — r^) + 2^' which displays 7 and A^ — T^ as (weakly) imaginary 
quantities. 

We still have to satisfy the two first conditions in (2.8). First of all, again using the 
Gauss constraint, we find that 
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Re{-A, sm{6) + F J = i^Re(-A,.E^ sm{6) + E'^T^) 
= ^Re{g + AyEy sin(5) - {E')' + E^T,) 

Ey G 

= — Rei^ + Ay smi6) + Ty) 

« ^Re{AySm{6) + ry) (2.10) 

where we have used the identity (E^)' — T^E^ ^TyEy = (which is due to the covariant 
Constance of the triads) and the reality of the triads. This demonstrates that the real 
part of the first quantity in the first line of (2.10) vanishes if and only if the last line 
in (2.10) is zero, at least on the Gauss-constraint surface. However, since Ax,Ay,6 are 
imaginary, the imaginary part of both expressions is already zero so that at least one of 
the fields A^., Ay, 6 would be expressible completely in terms of the triads and therefore 
the symplectic structure would become degenerate. This furnishes the proof. 
Since we are mainly interested in playing with a model that captures some of the 
features of general relativity rather than describing the cylindrically symmetric case,Q 
we take the following viewpoint : 

We want to preserve at least the reality of the diagonal part of the extrinsic curvature 
(the extrinsic curvature is diagonal for the cylindrically symmetric waves). Since the 
metric is manifestly diagonal, we can impose the simple reality condition that 

a-(3, A,-r,, A,, Ay (2.11) 



^The cylindrically symmetric case is not of physical importance anyway because an axisymmetric 
matter distribution is quite improbable 



are all purely imaginary. Now, looking at (2.8), we see that the extrinsic curvature will 
be symmetric upon imposition of the Gauss constraint but it will develop an imaginary 
off-diagonal part in the course of evolution. The only case when it is real is when the 
off-diagonal part vanishes and we are then back on a genuine subset of the phase space 
of the cylindrically symmetric model (this subset is a fix point under evolution). The 
diagonal part is real. The metric is diagonal and real by construction so we managed 
to impose reality conditions that are 'as close' to the cylindrically symmetric case as 
possible. 

These reality conditions are also preserved under evolution as we will show shortly and 
that raises another problem : Due to the evolution law of full general relativity 

qab = 2NKat + D^aNb) (2.12) 

an initially real, diagonal metric becomes complex and non diagonal when inserting 
our Kab above. This is apparently a contradiction ! However, there is a mistake in this 
argument : the processes of deriving the field equations and reducing down the degrees 
of freedom cannot be expected to commute if the phase space of the reduced model 
does not lie entirely in that of the full theory as it is the case for our model. 
One might ask whether it is possible to recover the cylindrically symmetric case by 
making use of all SO (3) degrees of freedom (the three Euler angles rather than only 
one, a and /3 for A^ and Ef respectively). However, if that was true one could always 
go to the gauge (2.2), (2.3) and since the reality structure of Einstein gravity is S0(3) 
covariant (since SO (3) is a real Lie group) one would obtain a contradiction. 
Thus, it is not possible to treat the cylindrically symmetric case while making use of a 
corresponding reduction on the Ashtekar phase space. 

Nevertheless, this model has all features in common with full general relativity (except 
for the modified reality conditions) so that it serves as an interesting new testing ground 
for the quantization programme of full quantum gravity as we will now show (actually 
this is also the viewpoint that one takes in dealing with 2+1 gravity as a toy model for 
3+1 gravity : the connection of 2+1 gravity is purely real). 

3 Complete symplectic reduction of the model 

It will turn out that the symplectic reduction programme can be carried out completely, 
that is, with respect to all constraints. 

After integrating over the (finite) range of the coordinates x and y, the reduced action 
becomes ([IH; we neglect a trivial prefactor coming from this integration) 

S= f dt f dz[-iA\Et - i^A'g^ - iN'^Va + iV^C)] (3.1) 

where it is understood that we have to substitute for Aj^, E^ the reduced expressions 
(2.2) and (2.3). Here Qi, Va, C are, respectively, the Gauss-, Vector- and scalar con- 
straint and A*, iV°, A^ are the corresponding Lagrange multipliers, called respectively 



the Gauss scalar potential, shift vector and A^ := Jdet{q)N the lapse function. 
By plugging the formulae (2.2) and (2.3) into the expressions for the constraint func- 
tions we obtain 

Gi := VaEt := daEt + e,,kAiE^, 



5,:,[{E^)' - sin(a - P){A,E- + AyEy)] =: 5,^g 



Vn := Fl.E, 



ah i 

= 6l[cos{a - P){{A^)'E'^ + {AyyEy) - A,{E')' - a'sin(a - P){{A^)' E'^ + {Ay)' E^] 
= 6:[cos{a - mA.yE^ + {AyYEy) - ^{E^y + a'G] =: <5^\/ 

/^ TTli 771a 7716 

= 2[A^AyE''Ey + {{A^yE'' + {AyyEy)E' sm{a - (3) 

+7(A^E^ + AyEy)E' cos{a - f3)] =: 1C (3.2) 

where F^^ is the curvature of the Ashtekar connection, ej^fe is the Levi-Civita symbol 
and we have abbreviated 7 := A^ + a'. Further we have absorbed the trivial factor of 
2 in the last line into the Lagrange multiplier of the scalar constraint. 
Now we insert the reduced form of our basic coordinates into the symplectic potential 
and obtain up to a total differential (note that spatial integrations by parts do not 
contribute boundary terms because the torus has no boundary) 

ie[(9i] = I dz[A,E' + A^E^ cos{a - (3) 

+AyEy cos(a - /9) - aiA^E"" + AyE^) sin(a - /5)] 
= / dz[{A, + a)E' + A^E"" cos(a - /?) + AyE^ cos(a - 13) + aQ] 

=: / dz[^ii^ + Aji^ + Ayiiy + dn"] . (3.3) 

Note that all configuration variables on the Gauss-reduced phase space are imaginary 
while the momenta are all real. This is a very nice reality structure. 
We next write the constraints in the so defined canonical coordinates. For the Gauss 
and vector constraint this is easy : 

Q = n^ 

V = {A^yw + {Ayyuy + a'w" - -f{wy (3.4) 

while for the scalar constraint we have to do some more work. 
We have directly from the definitions of our canonical coordinates 

tan(a - f3) = -^-^ — — (3.5) 

^ '^ A^U^ + AyUy ^ ' 

so that we have managed to express the combination a — j3 in terms of canonical 

variables. 

Next, we use the trigonometrical identities 



sm = ?r, cos 



tan^ ' 1 + tan^ 

and have unambiguously 

c = A,Ay—^^^^^^ + {{A^yT{^ + {Ayyuy)iri^n{a-^) + ^{Ajr + AyUy)iv 
= A^AyWWii + r (nj) - n" p. ((AJ'n" + (A„)'n^)m 



fn^)' - IT" 



Ajl"" + AyUy 



^ ''^A.^AyWW{[{Yr)' - n"]2 + [Ajr + A^^n^]^) 



+((Aj'n^ + (A,)'n^)n''[(n'')' - n"][A,n^ + A,n^] + 7(A,n^ + AyiiyfTV} 

We now rescale all configuration variables except for a by a factor of i so that these 
become real too (e.g. —iA^ will be called A^ again. Then the reduced action becomes 

S= f dt f dz[A^E'' + AyRy + ^YV + iaYV" - [iAH" + N'^V - NC]] (3.7) 

■J R -J Y^ 

with manifestly real constraint functions Q^ V and C and Lagrange multipliers. This 

furnishes the proof alluded to in the previous section that the reality structure of the 

model is preserved under evolution. 

We will assume that the prefactor 1/(^2:11^ + ^2^11^)^ never diverges or vanishes so that 

we can absorb it into the Lagrange multiplier and thus obtain a fourth order scalar 

constraint for this model. 

It is not trivial to check whether the constraint algebra closes. Before doing that, we 

will pass to the following equivalent set of constraints, obtained by subtracting a term 

proportional to the Gauss constraint from the scalar and vector constraint : 

Q = n", 

V = {A^)'ir + {Ay)'w-^{ir')' 

C = A.,AyU^W{-{{W)'f + {A.,ir + AyWf)+^{A^Il^ + AyIiyflV 

-{{A^)'W + {Ay)'W)W{ir)'[AJl^ + AyW] (3.8) 

and something amazing happens when adding the term n'''(n''')'[yl2:n^ + Ayll^jV pro- 
portional to the vector constraint to the scalar constraint 

c = A.,Ayii^iiy{-{{Yvyf + {Aji^ + AyYiyf) 

-^ir'[{IV)'f[A^Il^ + AyW] + ^{A^ir + AyWflV 
= [-{{ir)'f + {A^W + AyWf][A^AyWW + ^ir{AJl^ + AyW)] (3.9) 

i.e. the scalar constraint factorizes into two second order constraints with respect to 
the momenta ! They are given by 

Ci := A.^AyWW + -i A JY'W + -i AyWW audi (3.10) 

C2 := -{{Yr')'f + {A.M'' + AyWf. (3.11) 

By studying the transformations that these constraints generate, we obtain that 
Aa;, Aj^,7,7r^, n^, n''', n" are 0(2)-invariant while a is an 0(2)-angle and that 
Axi Ay, a, n'' are scalars while E^ , E^, 11", 7 are densities of weight one, i.e. Q generates 
0(2) transformations while V generates diffeomorphisms on the torus. By the way, 7 
can be interpreted as the abelian version of a connection plus the Maurer-Cartan form 
(^Mc = +{dS)S'^ (which in one dimension for S = exp(Q;) boils down to 9mc = Oi'dx) 



and therefore is an 0(2)-invariant quantity. 

The only nontrivial Poisson-bracket to check is then again that between two scalar 

constraints (all constraints are 0(2) invariant and Diff(S) covariant). We have 

{C[MIC[N]} = I dx j dyM{x)N{y){C,{x)C2{x),CMC2{y)} 



dx / dyM{x)N{y)[C,{x)C,{y){C2{x),C2{y)] + C,{x)C2{y){C2{x),C,{y)] 
+C2{x)CM{Ci{x),C2{y)} + C2{x)C2{y){C,{x),Ci{y)}] 

^jjxjjy{M{x)N{y)-M{y)N{x))[C,{x)CM{C2{x),C2{y)} 

+2Cr{x)CM{C2{x),CM} + C2{x)CM{Ci{x),C,{y)}] 
^dx J^dy{M{x)N{y) - M{y)N{x))Ci{x)C2{y) 

[-64x,y)]2{Wy{x)W{y){A,W + AyUy){y) 
= I dx{M'N - MN')CiC2i{Il^f)'{A^W + AyW) 
= C[{M'N - MN'){{Wyy{A^U'' + AyUy)] (3.12) 

i.e. the constraints form a first class algebra, open in the BRST-sense, and we can 
apply the framework of symplectic reduction (0]). 
We begin with the Gauss-constraint : 

We just have to pull back the symplectic potential to the 11" = surface, the vector 
and scalar constraint as they stand in (3.8) and (3.9) are already reduced. 
As experience with other model systems showsQ, it turns out to be much more conve- 
nient to first reduce the scalar constraint and then the vector constraint because the 
scalar constraint is not diffeomorphism invariant but only diffeomorphism covariant. 
We have to distinguish between two possible sectors : sector I is defined by the van- 
ishing of Ci, sector II by the vanishing of C2. 

We should mention that sector II corresponds to metrics that are degenerate in the 
z-direction : For our model the metric tensor takes the form 

^ EVE' E'E"" E'^Ey^ 
Qab = diag(-^, -^, -^) .. (3.13) 

From the Gauss-constraint and the definition of II^^^, IIj^ we infer 

{A^E^ + AyEyf = -[{Enr + [A.n^ + Ayliyr and |J = ^ (3.14) 

so that we can solve E^, Ey, E^ in terms of reduced coordinates 

(E", Ey, E') = (±-—-^^——-.fc'2, ±-—-^^——-\fc'2, m) (3.15) 

so that 

Oafe = diag(±— — — , ±— =^ — , ± , , ^ --I— rrr-— 62) . (3.16) 

^ 6V -Q^. ^y [A.^Yl'^ + AyW^W ' ^ ' 

Note that the quotient II^/II^ remains finite. 
This shows that sector II is the 'unphysical' one. 

^specifically, the author was looking at spherically symmetric gravity plus matter (P|) 



3.1 Sector I 

We pass to new canonical pairs, 

(g. := H\A,\),p^ := A^W; g, := ln(|A,|),p^ := A,n^; q := ln(|7|),p := jU^) (3.17) 



(note that 51n(|a;|) = 51n(a;) for real x where |a;| := vx^) because then the scalar 
constraint Ci adopts the simple form 

p^pv ^ pVp ^ pp^ = Q (3.18) 

which does not even make any ordering problems any more if one would pass to the 

quantum theory directly, that is using the Dirac approach. Moreover, it defines trivially 

an abelian constraint because it consists only of momenta. 

Note that we could have changed the polarization for the pair 7, 11'' in order that 

all configuration variables are scalars without destroying the simple form of the scalar 

constraint, but then we would leave the Ashtekar-polarization. Let us therefore proceed 

this way although this will equip the canonical coordinates with a somewhat awkward 

tensor valence. 

It is clear that one will now diagonalize this constraint because then one can solve the 

Hamilton- Jacobi equation by separation of variables. We have 

C = ^[(p" +P^ + 2p)2 - (p" -p^) - (2p)2] (3.19) 

so that we define still another set of canonical pairs 

P° := lip' + p' + 2p), P' := \{p' - py), P' := p (3.20) 

and conversely 

px ^pO^pl _ 2p2^ pV = pO_pl_ 2p2^ p ._ p2 (3 21) 

which allows for the determination of the new canonical configuration variables via 
inserting (3.21) into the symplectic potential 

e[dt] = |rfa:fe(P° + P'-2p2)+g,(P°-pi-2p2)+gp2] 

= Jjx[P^j^{q, + qy)+P'j^{q,-qy) + P'j^iq-2q,-2qy)] (3.22) 
from which we read off 

Qo = qx + Qy, Qi = Qx- Qy, Q2 = q - 2q^ - 2qy . (3.23) 

Then the scalar constraint just says that (P°, P^, P^) is a null- vector in a 3-dimensional 
Minkowski space at every point of S. There are two branches of general solutions to 
the associated Hamilton- Jacobi equation, 

S±[p'',p'';Q,] = [ dz[±^{p-Y + {p^YQ,+p''Qi+p'Q2] (3.24) 

10 



where, following the second reference of ^, Pu = P^,Pv = P^ are the new momenta 
and 

q. := ^ = ± , ^^ Q' + Q' (3.25) 

are the new invariant (under the motions generated by the scalar constraint) configu- 
ration variables. 

One can get them also just by pulling back the symplectic structure to the scalar con- 
straint surface via its embedding l into the phase space (compare 0). Inserting the 
solutions P° = ±J{P^y + {P'^Y iiito the symplectic potential yields (as before we do 



not display total differentials) 

= - f r(± ^" ^Qo + Qi) + p"(± , ^' =go + Q2) 

= I b"^(± , ^" Q, + Q,)+p^^{± ^^ go + ^2) (.3.26) 



It remains to reduce the vector constraint. Resubstituting all the symplectomorphisms 
we have carried out up to now, we obtain (note that for any two functions a and b 
holds 



'a+{ , ^ :yb = ' , ' ^J ^- ; , ' i_ (a^ + b^) = 



Wa^ + b^ Wa^ + b^ 2Va^ + b^ 2(V^^TF)= 

V = {A^yw + {AyYuy - j{Ptry 

= fe)y + (g.)y-7(— )' 

7 

= (q^yp^ + (qvYp'^ + q'p - p' 

= {Qoyp'+{Qiyp' + {Q2yp'-{p'y 



±(go)V(p")' + ip")' + iQi)y + (Q2)y - iPi' 
{quYp'' + (g.)y - {p^y 

(gjy + (ln(|g^|))'g^p--(g^p-)' 

(g«)y - qn,{pn' (3.27) 



where we have carried out one more canonical transformation 

q, =: ln(|g^|), p" =: q^p"" . (3.28) 

Equation (3.28) displays qu,p^ as scalars and qw,P^ as densities of weight one. 
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3.1.1 Reality structure of sector I 

From their definition (3.20) it is clear that p^,p^,p as products of two real quantities 
are all real, whence also P^ are real. 
From (3.17) we have that 

go = H\A.A,\), Q, = H\^\), Q, = ln(|— 1^1) (3.29) 

SO also the Q^ are real. Therefore, altogether, the left over canonical pairs {qu,P^', (lv,p") 
are both real and have range over the whole real axis, whereas due to g^ := exp(gt,) 
has only positive range while p"" := p" /qv has still range over the whole real axis. 

3.2 Sector II 

Now the scalar constraint C2 can even be written as one of two branches of a constraint 
linear in momentum 

C± := Ajr + AyW ± (n^)' . (3.30) 

Quite similar as for sector 1 we define new canonical pairs 

(g, := H\A^\),p^ := A.,W- g, := ln(|A,|),p^ := AyW- q := -f,p := W) (3.31) 

and obtain 

C±:=p''+py±p' . (3.32) 

The general solution of the Hamilton- Jacobi equation is obviously given by 

Si = Jjx[{p^yq, + {p-yqyT{p''+p-)q] = jjx{p\Tq-{q.)')+p'"{Tq-{qy)')] (3.33) 

whence we read off for the new momenta 

{Pl' = P\ ipn' = P' (3.34) 

and for the new configuration variables 

qv ■= Tq - {qx)\ qw ■= Tq - (qy)' (3.35) 

which are both densities of weight one so that the reduced (with respect to the scalar 
constraint) vector constraint becomes 

V = -q.ipn' - q^py = iqu)y - q^ipy (3.36) 

where we have defined still another canonical pair by 

g„ =: (g„)', g« := -{p^)' . (3.37) 

That this defines indeed a symplectomorphism is obvious from the following short 
calculation (spatial surface integrals vanish due to our choice of the spatial topology) 

9.p^= /a)y= f u-pn'- (3.38) 



12 



3.2.1 Reality structure of sector II 

Again we liave tliat p^,py,p are all real which then is true via (3.34) also for p^^p" and 
via (3.31), (3.35) also for qv,qw Finally, via (3.37) also for qu,qw This time, however, 
g^ is not necessarily positive. 



3.3 Reduction of the vector constraint 

We come now to the main problem for a spatially diffeomorphism invariant field the- 
ory, namely to find a complete but minimal set of Dirac observables which are to be 
promoted to the basic set of quantum operators later. 

The stress of the following exposition is on the conceptual side. It is not meant to be 
complete or rigorous to the last detail. 

It is always easy to find (classical) expressions that Poisson-commute (even strongly) 
with the diffeomorphism constraint - the recipe is the following : 

1) From the transformation law of the basic field variables under diffeomorphisms de- 
rive their tensor nature. 

2) Construct scalar densities of weight one. 

3) Integrate them over the initial data hypersurface. 

These objects are then already invariant under diffeomorphisms that are connected to 
the identity (in case of an asymptotically fiat topology in general only if the diffeomor- 
phism tends to the identity rapidly enough at spatial infinity). This is obvious from 
the transformation property of a scalar density s of weight one under an infinitesimal 
diffeomorphism (generated by a vector field C,"') 6s = da{C,"'s). 

However, the scalar constraint is not at all so easy to solve since it can be viewed as 
a dynamical constraint. Hence, it is most important to keep the algebraic structure of 
the scalar constraint as simple as possible. 

Now, since the scalar constraint is a density of weight two, it has a nontrivial transfor- 
mation law under diffeomorphisms. Accordingly, solving the vector constraint before 
solving the scalar constraint will most likely complicate rather than simplify the alge- 
braic structure of the scalar constraint because it will become a non-local object and, 
moreover, cannot be reduced because it cannot be written in terms of diffeomorphism 
invariant quantities. 

The conclusion of all this is that the scalar constraint is to be solved before the dif- 
feomorphism constraint for a field theory with a scalar constraint which is at least 
quadratic in the momenta. This is what we were able to do up to now. 

Now we are going to derive a sequence of results which generalize in an obvious manner 
(although this will become combinatorically much more difficult) to higher dimensions. 

Lemma 3.1 The set of all integrated scalar densities constructed from the configu- 
ration variables provides for a (possibly overcomplete) system of coordinates of the 
diffeomorphism reduced configuration space. 

Proof : 

1) Any local object constructed from the configuration variables has a nontrivial trans- 
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formation law with respect to infinitesimal diffeomorphisms. Accordingly, diffeomorphism- 
invariant objects are non-local, i.e. they can be expressed as integrals over (parts of) 
S of local objects. 

2) Were these integrals not over all of S then they were not invariant because the 
diffeomorphisms have their support everywhere except for the boundary of E. 

3) Only scalar densities of weight one are invariant when integrated over S. 

4) One could consider objects of the form 

/ dxi... / dxmK{xi,...,Xjn) (3.39) 

i.e. K[xi, ..,Xm,) is some integral kernel. However, any such integral kernel can be 
approximated arbitrarily well (in a given suitable norm by an argument along the lines 
of the nuclear theorem) by linear combinations of tensor products fi{xi)..frn{xm) (for 
each k the different functions fk will belong to a base). Thus, the integral of K can be 
written as a sum of products of single integrals which displays K as a derived object. 
This proves the lemma. 
D 

Let {q^}i(zi denote this set of Dirac observables where the labelling is with respect to 
an index set I whose cardinality is to be obtained. Let us further assume that all these 
observables are algebraically independent. Any linear combination 

S■.= J2P^^' (3-40) 

of these observables is again an observable where the pi are, in general, complex num- 
bers (the reality conditions on the pi match those on the g* so that there are as many 
independent momenta as configuration variables). Then, provided that the set of the g* 
is complete in the sense that every diffeomorphism invariant object constructed from 
the configuration variables can be expressed by them, the above functional S is the 
Hamilton- Jacobi functional relative to the diffeomorphism constraint and the Pi play 
the role of the integration momenta (recall the formalism in the second paper of 0). 
We will now define an overcomplete set of diffeomorphism invariant functionals of the 
configuration variables for our model and then systematically obtain a minimal subset. 
When stripping off the technical methods that are particular for this model, one will 
get an insight for what is necessary for the full theory. 

According to lemma (3.1) what we have to do first is to construct all possible scalar 
densities s that can be built from Qu, qw We first need two technical 
lemmas : 

Lemma 3.2 1) In one dimension, tensors of valence (r,s) are scalar densities of weight 
s-r. 

2) In one dimension, a scalar density s of weight one gives rise to a s-compatible 
connection T according to 

r := [ln(S)]' . (3.41) 

Proof : 

1) 

This lemma is most easily proved by comparing the corresponding transformation 

properties of tensors T°-' b.. under diffeomorphisms x" — > x" = x" — ^" 

5T"- b.. = L^T'^- ,.. = CT"- b..,c - T^- b..Cc - •• + T''- c.Cb + - (3.42) 
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which in one dimension reduces to 

5T = iT' + (s - r)i'T . (3.43) 

On the other hand, a scalar density s of weight n transforms as 

5~s = rs,a + n~sil (3.44) 

which again in one dimension reduces to 

5~s = ^§' + ns^' (3.45) 

and that furnishes the proof of the first part of the lemma. 

Accordingly, in one dimension tensor valences are just characterized by their density 

weight. 

2) 

One simply builds the unique torsion-free connection compatible with a given metric 
by restricting the Christoffel formula F^^ = l/2q"''^{qdb,c + (ldc,b — (lbc,d) to one dimension. 
A metric in one dimension is, according to 1), just an arbitrary density of weight 2 
corresponding to a symmetric tensor of valence (0,2) (i.e. a metric q). Since the inverse 
of a metric in one dimension is just one over this metric we have for the afiine connection 
in one dimension F = l/2(g + q — q )/q = l/2[ln(g)]' or, by defining the density s of 
weight one through the equation q = s^, we have finally a connection defined by an 
arbitrary density s of weight one given by 

F := [ln(5)]' . (3.46) 

Again, by just restricting the covariant derivative 

D,T«- b.. = d,T^- b.. + Tl^T"- b.. + .. - Ff,T"- ,.. - .. (3.47) 

to one dimension we find that for a density T of weight n we have 

DT = T' - riTT (3.48) 

and one can explicitly check that DT is a density of weight n+1, i.e. a density again 

with one more covariant valence. 

In particular Ds = 5' — Fs = which shows that the connection is s-compatible. 

D 

We are now ready to prove one of the main theorems of this section. 

Theorem 3.1 An (overcomplete) system of scalar densities of weight one is given by 

9-/({(g«)^"^}neA^) (3.49) 

where f is an arbitrary function of the infinite number of arguments indicated, M is the 
set of non-negative integers and for any scalar function s, s*^") denotes the scalar d^^s) 
where d is the operator d = l/q^d which is obviously anti- self- adjoint with respect to 
the measure djj, = q^dz. 
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Proof : 

It is clear that any scalar density s of weight one can be written as q^s where s is a 
scalar, simply define s := s/qw. The proof thus boils down to showing that the above 
function f in (3.49) is the most general scalar that we can build. 

Such a scalar will be the most general scalar constructed from g^, g^ and all their par- 
tial derivatives up to arbitrarily high order. This can be done in a covariant way only 
by contraction of tensors, that is in one dimension, by multiplying densities such that 
their overall weight is zero, according to lemma (3.2), first part. Hence we need to 
construct a covariant derivative such that the derived tensors, i.e. densities, transform 
as tensors, i.e. densities again. This can be done by constructing an affine connection 
in one dimension and here part 2) of lemma (3.2) comes into play. 
We are now able to take arbitrarily high covariant derivatives of qwiQui for example 
with respect to the connection derived from s := g^. But then Dq^ = because D is 
Qw compatible which shows that the covariant derivatives of q^ do not contribute in our 
construction and furthermore D"-(f) = g^(l/g^D)"0 = ql^d"'^) where we have introduced 
the operator d = 1/qwd. 

It remains to show that there is no loss in generality in choosing s := g^. 
Assume that we had chosen any other density s instead of g^ in the construction 
of r then the covariant derivatives D of g^ would not vanish, however then Dq^ = 
sd{qw/s) = sD{qw/s) = —qw/sDs and s is again built only from our configuration 
variables and D (it must be a tensor and it cannot be constructed by using D because 
D is defined through s). Hence Dq^j can be recast in the form of a function f of the 
type given in (3.49). 
This furnishes the proof. 
D 

What we have achieved up to now is to obtain a sufficient number of observables g* 
which are functionals of the configuration variables only. 

We now ask for possible redundancies among them. Two sources of redundancy are 
obvious : 

1) by doing integrations by parts it is possible to relate some of the g* and 

2) there might be topological identities which could arise, for example, via index the- 
orems and thus affect the range of g* (for example in 2 dimensions the integrated 
densitized curvature scalar is an object of the form we are considering here, but the 
Gauss-Bonnet theorem states that it depends only on the genus of the compact hyper- 
surface without boundary). 

Under the assumption that we can restrict to analytical f 's in (3.49) it is possible to get 
rid of the ffrst redundancy. What we are presupposing is that there exists an expansion 
of f of the type 

/({omlmeA^) = Y. /{"fc} n Ofc' (3-50) 

{{nfc}} fe=0 

where a^ := d^qu-, the sum is over all sequences of (non-negative) occupation numbers 
Uk and the coefficients of the monomials 

oo 

M{n,} ■■= n «r (3.51) 

fc=0 

are real numbers because g„, q^ are real as shown above. Hence, in our case the index 
set is given by the set of sequences, I = {{n^}} 3 i, and the (still overcomplete) set of 
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g*'s is given by the integrals 

q' := I dzq^M' (3.52) 

whereas the fi are the conjugate momenta. 

We will divide out the first, combinatorial redundancy by bringing the g* into a stan- 
dard form thereby introducing a terminology that one is used to from Fock-space or 
statistical mechanics. : 

First we observe that by doing an integration by parts the occupation numbers n^ of 
the quantum states or energy levels k for particles of type g„ may change but that the 
total particle number (of each particle type separately if there are more than 2 fields 
on the configuration space which it is not the case in our situation) 

oo 

N:=Y.n, (3.53) 

fe=0 

as well as the total energy 

oo 

E:=Y^ krik (3.54) 

are preserved under such a process, simply because the number of factors of a (irre- 
spective of how often they are derived) and the number of derivations does not change. 
In fact, one could write the operation of doing an integration by parts in terms of 
annihilators and creators. 



Lemma 3.3 Consider now the class of (g*) 's with the same total particle number N 
and total energy E. Then an algebraically independent genuine subset of g* 's in this 
class is given by those with indices i such that rriQ = l..N,mi = 0, the rest of the 
occupation numbers being arbitrary (subject to the condition that they lie in the class 
labelled by (N,E). 

Proof : 

In the first step we show that every member in the given class can be written as a 
linear combination of these g"s and in a second step we show that these g* cannot be 
transformed into each other. 

Consider the class (A^, E) = (M, k) and we may exclude the trivial case M = 0. Then 
g* is of the form 

g' = / dfi aki..aE^(ki+..+kM-i) (3.55) 

where we do not care about occupation numbers, i.e. some of the kr may be equal. Let, 
without loss of generality, ki be the lowest of the kr. If /ci = we are done, otherwise do 
an integration by parts. After that, g* is a linear combination (with integer coefficients) 
of g-^'s of the required form since all these q^ have mo = 1. By doing further integrations 
by parts if necessary for these q^ (since k{ao)''~^ai = (i((ao)'^) can be integrated by 
parts) one can satisfy mi = 0,mQ > 0. This finishes the first step. 
Let now one of these 'reduced' g*'s be given. Then the set of its occupation numbers 
is characterizing, i.e. by doing an integration by parts, starting from mi = 0, one 
necessarily picks up terms that have mi ^ from the differentiation of a since we 
already have tuq > 0. This furnishes the proof. 
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n 

From now on we will understand the index set I to be reduced, that is, we consider 
the 'Mandelstam-identities' between the original g*'s originating from integrations by 
parts to be taken care of by lemma 3.3 . 

As far as the second source of redundancy is concerned, we do not have any indication 
that any of the g* should have a non-continuous range. 
There is still a more subtle source of redundancy : 

One could fix a gauge (a frame) and expand g„ into a taylor series with respect to 
this frame where we have assumed these variables to be smooth functions on E. The 
freedom in choosing these functions is then captured in the choice of the values of their 
taylor coefficients. 

However, then one could just integrate out all of our reduced g* with respect to this 
frame and we would get complicated rational functions of these Taylor coefficients 
multiplied by real numbers coming from the integration of the powers of the frame- 
variable X over the hypersurface and depending only on the topology of S. 
Since any such rational function can be generated from the Taylor coefficients itself, the 
knowledge of the latter is sufficient and provides for the minimal subset of observables. 
The result of all this analysis summarized in the subsequent theorem. 

Theorem 3.2 Only the following g* 's are independent variables on the dijfeomorphism- 
reduced configuration space (n is any non-negative integer), possibly modulo some dis- 
crete degrees of freedom : 

qn--=l f d^ud^'^qu = (-1)"^ / rf/i[rf"g„]' . (3.56) 

Proof : 

There is a one to one mapping between these functionals and the Taylor coefficients of 

g^, possibly up to sign. For example, after fixing a gauge, e.g. g^(a;) = g'{x) = 1, d = d 

where g is a fixed scalar function, the g„ allow to extract all the Taylor coefficients of 

Qu (e.g. d'^cj) starts with the n-th coefficient), possibly up to sign (that is, we might 

miss some discrete degrees of freedom) . 

D 

Note that (3.56) are jA/"! independent quantities (A/" = natural numbers), where \S\ 

denotes the cardinality of the set S, because they satisfy all the criteria of the previous 

theorem to be algebraically independent and precisely for that reason we have 6„/2 = 

$n/2 = 0. Accordingly, the number of degrees of freedom is naturally countable without 

employing a complete system of mode functions on S to express the true degrees of 

freedom in discrete form. 

Moreover, these observables have the nice feature that 



qn,d'''qu{x), (3.57) 



j2n 

, - = qwf- 



-%T = -\^ {-l)\d\u{x)){d^--\u{x)) . (3.58) 



We will denote the momenta conjugate to g„ by p^' 
The associated Hamilton- Jacobi functional is 



-5:=^^"^" (3-59) 



n=0 
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and the old momenta are the functional derivatives of S with respect to the old coor- 
dinates. 

For sector I we have seen that g^ is positive (and therefore serves as a genuine density 
for a measure on S in the the definition of Qn) whereas for sector II the range of q^ is 
the whole real line. Since g„, g^ are both real, {d^QuY is positive whence g„ has half 
or the full range of the real line for sector I or II respectively, p" on the other hand 
is generally real so that the pair (g„,p"') coordinatizes the cotangent space over the 
(positive) real line. We can avoid this complication for sector I by going to the pair 
{Qn,P^) -^ (}^{\Qn\) , InPn) which we will assume in the sequel. 

Our, admittedly somewhat sketchy construction generalizes in an obvious way to more 
dimensions for every difFeomorphism invariant field theory. 

4 Quantum theory and determination of the phys- 
ical inner product 

4.1 Reduced phase space approach 

Quantum theory is straightforward. Keeping the polarization that we have ended up 
with, physical states depend only on the variables g„. 

It is important to see that the scalar constraint forced us to let the reduced configu- 
ration space (with respect to the scalar constraint) not only depend on the original 
configuration variables (before reducing) but also on the original momenta ! This mix- 
ing of original polarizations is typical for genuinely quadratic constraints with respect 
to the momenta and can be viewed as a source of difficulty in the Dirac approach to 
select physical states. 

The Dirac observables qn,P^ are both real and have both infinite range according 
to the agreement that we have made above. The unique inner product that turns 
qn,p"' into self-adjoint operators (for a specific choice for n) would be the one corre- 
sponding to the Lebesgue measure dqn- Since we are dealing with functions that are 
defined on the infinite direct product space x'^^qR we have the problem that the for- 
mal "infinite product Lebesgue measure" A^o'^'?" which would turn all the qn,p"' into 
self-adjoint operators is ill-defined unless we integrate only cylindrical functions. A 
solution is to go to the Bargman-Fock representation and to work with the observables 
o-n '■= Qn + ip^ and a„. As is well known, the measure /i that implements the reality 
conditions (a„)^ = a„ is Gaussian and therefore we can actually integrate holomorphic 
functions on X := x^qC that depend on an infinite number of variables. We therefore 
end up with L2(X, /x), the usual Bargman-Fock Hilbert space. 

4.2 Dynamics of the model - deparametrization for an infinite 
number of degrees of freedom 

What we have obtained up to now is the complete set of Dirac observables which are 
time-independent. However, their interpretation is still missing. On the other hand, 
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interpretation is usually straightforward for non-gauge invariant quantities. In order 
to arrive at an interpretation there are two possible strategies available : either one 
expresses an object that is 0(2) x Dif fill) invariant in terms of Dirac observables 
plus an additional gauge function t which labels the choice of foliation of spacetime 
and quantizes only the Dirac observables (|ll|) or one uses the framework of time- 
dependent Dirac observables (deparametrization, |^). 
Let us try to apply the latter one. 

In the sequel we will only deal with the 'physically relevant' sector I. 
We are asked to cast the scalar constraint into a form that allows to apply an extension 
to an infinite number of degrees of freedom of the framework of deparametrization. As 
is proved in ([Q) the only change is that the time variable becomes 'many- fingered' 
(Tomonaga-Schwinger extension). 

In particular we are interested in the time-development of the quantities E^,Ey,E^ 
which determine the metric. 

The starting point will be equation (3.19). We take the square root and restrict to the 
positive branch. Then the scalar constraint adopts the familiar form 



P'{z) - H{Q,, P'; Q2, P'){z) := P%z) - ^[P\z)f + [P2(^)]2 (4.1) 

which is already deparametrized, although at each point x of the hypersurface sepa- 
rately. In (4.2) we have defined an infinite number of 'Hamiltonians', H{x), that is, 
an 'energy density'. Note that, since H involves only momenta, there are no ordering 
problems. 
Following this framework, we construct the time-evolution operator 

Ut := exp(-« I dziQo - t)H) (4.2) 

and define any operator O (to be factor-ordered appropriately) on the set of physical 

states ^ = Uip (with respect to the scalar constraint) by 

{dt^)[Qo,Qi,Q2] := UA=qMQ\Q'] ■ (4.3) 

Now, in contrast to toy models for cosmology, physical states should also be diffeomor- 
phism invariant. We will therefore take the operator to be diffeomorphism-invariant, 
but not necessarily invariant under the scalar constraint. 

We want to apply the above framework for example to the operator for the three- volume 
whose classical counterpart is given by 

V := f dz^Jdet q = f dz\/WEyE^ . (4.4) 

We will proceed as follows : We start from (4.2) and proceed to the Dirac quantization 
of the model. As for the scalar constraint, we have no operator-ordering problems and 
the complete solution of 

(_[p0]2 ^ [pip ^ ^p^f)^,yQ^^Q^^Q^] = (4.5) 

is formally given by a functional integral (positive energy solutions) 
^Qo,QuQ2] = J ^ [d^k](f)[k\ P] exp(-z J dz[^ {k^f + {k^YQo - k^Qi - fc'QJ) 
= exp(-2 J^ dzQo^iP^r + (P^rmQi, Q2] (4.6) 
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where [d'^k](f)[k^,k'^] is the formal expression for the hmes of a cyhndrical measure 
(formally, [d'^k] is an infinite Lebesgue measure and is the Fourier transform of \E'). 
The integral is defined through a lattice evaluation and taking the limes that the lattice 
spacing goes to zero at the end (the lattice is finite since S is compact). 
Since Ut commutes with the vector constraint (ordered with the momenta to the right), 
diffeomorphism invariance is obeyed if and only if for an infinitesimal density C, of weight 
-1 holds 

mi + (Qi)'e, Q2 + (^2)'^ + e] = i'lQi, Q2] ■ (4.7) 

Inserting, the integrand of the path integral becomes after doing an integration by 
parts in the exponential (up to first order in ^) 

0[A;\ e] expit f dx[ik' - {^k'y)Qi + (k' - i^k^y)Q2 + Ck']) ■ 

We change to new variables K^ := /c* — (^A;*)', i = 1,2 in the functional integral and 
obtain (up to first order in ^) the condition 

(P[K, + {K^y, K^ + {K'O'] exp(^ / dz^K^) = m\ ^'] (4-8) 

where we have used the diffeomorphism invariance of the Lebesgue measure (the Jaco- 
bian is 1 to first order in ^ : denote by h the spacing of the lattice, k = k^ or k"^, kn '■ = 
k{nh), then the Jacobian is 

dkn d[K^ - l/{2hMK)n+l - (e^)n-l)] 



dKm OKra 

^n,m TTTTv ™i"*~l "n,m+lj • •Jnm • \^ ~r ^{C,} )n,m v^'^^J 

where the matrix A is trace-free and of first order in ^, so det(J) = 1 + ii{A) + .. = 

i+o(a)- 

We now split (t)[k^, /c^] =: expii ^-^ dz f {k'^))ip[k^ , k"^] and seek to determine the function 
f according to the requirement that 



J^ dx[f{K' + i^K^y) + ^'K'] = J^ dxf{K' 



up to first order in ^. We expand f and perform integrations by parts until only ^ 
appears derivated with respect to z. Then a sufficient condition for the last equation 
to hold is that 

The solution of this ordinary first order differential equation is given by 



f{k') = [c-\n{k')]k' 

where c is a real integration constant. 

This means that il) is diffeomorphism invariant when A;^,A;^ transform as densities of 

weight one. Again we are therefore lead to the functionals 

^[A;i, k^] = /({g"}), where g" = / dii^d'^^k^ /k^)f (4.10) 
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with dfi = dzk^, d = 1/k^d. 

Diffeomorphism invariance up to first order implies diffeomorphism invariance with 

respect to the component of the identity of the diffeomorphism group. 

We are now in the position to define the operator V. 

By making use of the definitions we can express (4.5) in terms of Q^, P'^ 






det(g) = exp(— Q 



0^ 



(pO _ p2)2 _ (pl^ 
2(P0 _ p2) 



4(po _ p2)2 _ [(exp(-(Q2 + 2Qo))P^)f (4.11) 

and we face severe factor ordering problems when promoting the integral over this 
function to quantum theory. 

The action of this operator on physical states is given by (we assume that we order all 
momenta to the right hand side in the formal power expansions of the non-analytical 
terms in (4.12)) 

VMQo, Qu Q2] = tJt J[d^k] exp(z J dz[Qik' + Q2k^ +[c- \n{e)]k^])i>[k\ A;^] x 

2 _ 



r /(A;2)2 _ (A;l)2 

j^dxexpi-t y'l ' ^4{k')2-[(exp{-iQ, + 2tW)f 



(4.12) 

which, when acting with the evolution operator (which, by the way, is unitary for 
Qo — t =const. with respect to the inner product (4.1) : since the integrated Hamilto- 
nian commutes with both constraints, it is physical, that is, well-defined on Tiphysi ^-nd 
since it is classically real it promotes to a self-adjoint operator by a suitable choice of 
ordering) results in a horrible object which we cannot work out explicitly. 
The point of this subsection was to show how deparametrization works in principle for 
field theories and how one can arrive at an interpretation. It also shows that most of 
the evolving operators which one would intuitively think of as observables are probably 
not well-defined on the physical Hilbert space (here in the connection representation). 



5 Loop variables 

The concluding section of this paper is designated to the issue of using (traces of) 
Wilson-loop observables for the quantization of or even for the classical treatment of 
gravity. Because our model system is completely integrable one can ask and answer all 
kinds of questions that have been raised in the literature. 

In particular we ask if the operators that are usually defined for the loop representation 
in the hterature ([^) are well defined at all, that is, if they make sense when applied 
to physical states. Of course, as they stand they will not leave the physical Hilbert 
space invariant and thus are not even expected to have a well-defined action on physical 
states. However, it will be interesting to see how the Dirac observables look in terms 
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of loop-coordinates, in particular the constraints, what their Poisson structure is and 
how one could define a 'loop-transform' into the loop representation (0), from states 
that have been defined for the connection representation. 

Also, one could construct loop-operators as certain 'time-dependent' Dirac observables. 
We will deal here only with the definition of appropriate loop variables in the connec- 
tion representation. 

We begin by recalling that one defines the following 'basic' loop variables when dealing 
with the loop representation (starting from the connection representation) : 
The T° variables are defined to be traces of the holonomies (parallel transport opera- 
tors) 

U^[A]{x):=Vexp{£A) (5.1) 

around closed curves (loops a) with base point x in the fundamental representation of 
SU(2) (the Ashtekar variables are complex-valued whence one they lie in the Lie algebra 
of complexified SU(2,C)=SL(2,C); although the gauge group of gravity is really SU(2), 
the loop variables are even invariant under SL(2,C)) 

T^A] := hriUM]] (5-2) 

where V denotes path-ordering and A := A^rfx^Tj, Tj a basis in the Lie algebra su(2). 
We will choose Tj := —i/2ai, the latter being the usual Pauli matrices. Note that the 
trace of the holonomy is independent of the base-point on the given loop. 
The T^-variables are defined to be 

T:[A,E]{x) := tT[E%x)UM]ix)] ■ (5.3) 

As one can show, they form a closed Poisson-algebra and can thus serve as a starting 

point for the group-theoretical quantization scheme. Let us explore what happens with 

these variables for our model. 

The (pull-back under a section of the principal SU(2)-bundle of the) connection reduces 

to 

A{x, y, z) = A^(cos(Q;)ri + sm{a)T2)dx + Ay{— sin(a)ri + cos{a)T2)dy + A^r^dz (5.4) 

and the gauge transformations generated by the Gauss-constraint for our reduced con- 
figuration variables are those under the U(l) subgroup of SU(2) generated by t^, (even 
its complexification) : 

- dSS-^ + SAS-^ where S = exp(Ar3) (5.5) 

can be written as in (5.4) only that a, A^ are replaced hj a + A, A^ — A'. Under a 
gauge transformation, the holonomy and the triad E = E^Tida transform covariantly 

E{x) -^ S{x)E{x)S-\x) and f/^(^) -^ S{x)U^{x)S-\x) . (5.6) 

whence the T^,T^ variables live on the Gauss-reduced phase space. 

It is convenient to eliminate the angle a by doing an U(1,C) rotation with gauge 
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parameter A = —a (note that a is complex so this transformation is really an element 
of U(1,C)). Then the connection and the triad adopt the reduced structure 

A{x,y,z) = A^{z)ridx + AyT2{z)dy + -f{z)T3dz (5.7) 

E{x,y,z) = [W{z)d,-Ey{z)sm{a-P)dy]n + [W{z)dy-E-{z)sm{a-(3)d^]T2 
+W{z)Tsd, . (5.8) 

Note that both quantities do not have any x,y dependence ! 

We will now start to find suitable loops to capture the full U(l) invariant information 
needed in order to actually coordinatize the Gauss-reduced phase space. Three loops 
are obvious : they are those lying entirely in one of the three 5'^-factors of the three- 
torus. We will call them a{x),a{y) and a{z) respectively. A short calculation yields 
for the corresponding T°, T^'s (the range of the three coordinates x,y and z is taken to 
be the interval [0, 1] with endpoints identified) : 

Ta{x){z) = cos{-A^{z)) (5.9) 



Ta{y){z) = C0S{- Ay{z)) 
1 /•! 

T^(^) = cos(-y dz-f{z)) 

T:i.){z) = -U-{z)smi^A^{z)) 

7^4) W = -Tl'{z)sm{^-Ay{z)) 

T:^,^{z) = -Uy{z)sm{y\z^iz)) 

T- (z) - ("-"')(^)° ,in(\i(z)) 

rr.v / X (n„nT)(2) ,1 , ,, 

^^(^)(^) - iA.l\AyW)iz) ''<-2^^^'^^ 

t:^') = Ti^^z) = t:^.4z) . = T:^y)iz) = o 



5.10) 
5.11) 
5.12) 
5.13) 
5.14) 
5.15) 

5.16) 
5.17) 



Note that from the four quantities To-(^), T^/ Jz), fj, G {x, y} we can already extract the 
two gauge invariant canonical pairs A^, 11^ locally while the two variables T^/^-s{z),Tcr(z) 
allow to recover 11''' locally but only global information about 7 through the integrated 
quantity Jj^dz'y{z). This has two consequences : 

1) The information contained in the two 'off-diagonal' variables Tj^t^Jz), fi ^ u is 
already redundant and so we can neglect them, 

2) We need a further T^ variable in order to get hand on 'y{z). 

Before doing that we show that the variables defined above form a closed (classical) 
*-Poisson algebra (the Loop- variables for the Ashtekar phase space are not closed under 
complex conjugation at least it is not obvious how to express T^ = |tr('Pexp(/^[2r — 
A])) in terms of T^ ^T^ where F is the spin-connection). 

We have proved in the sections before that A^., Ay, 7 are purely imaginary while 11^, 11^, 11^ 
are purely real. This implies the *-relations 



T.(„) (z) = T„^a) {z) and T;(,) {z) = -T^^,^ {z) (5. 1^ 
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i.e. the T°'s are purely real, the T^'s purely imaginary, and also that the range of 
the T° are the real numbers larger than or equal to one while the T^'s have the whole 
imaginary axis as their range. This opens the chance for a cotangent topology of the 
associated phase space and therefore a group theoretical quantization scheme. 
Another consequence is that the variables To-(^) are not quite sufficient to recover the 
connection up to gauge transformations : we have the degeneracy 

T.(^)[A]=T.(^)[-A] (5.19) 

which means that these loop variables are not separating on the moduli space A/Q 
of SL(2,C) connections modulo gauge transformations for our model when restricting 
ourselves' to the above defined class of loops. This is in analogy to the known degen- 
eracies ([|13|) of these loop variables as coordinates for A/Q^. 

Next we compute the Poisson-structure which turns out to be quite similar to the one 
obtained in the first reference of (|P) when restricting ourselves to the 
z-direction : 

(5.20) 



{Taia)iz),T^(b)iz')} = 


= 


{t.(,)(^),t;(,)(/)} = 


= ^5,Jiz,z')sin'iU,iz))- 


{T.az),T:^4z')} - 


. ls.n^lj\iz)dz)^lil 


R"(,)W,^;h(^')} = 


- {T.'a^),T:^.)i^')} = 


{t.%)(^),t;(,)(/)} = 


= 2 [^<^%) ^^'> " ^kz) (^')]^^(^) 



(5.21) 

[T.i.)?) (5.22) 

(5.23) 

(5.24) 

and we could even have written the right hand side of all the equations as linear 
combinations of single T°'s or ,T^'s with distribution-valued coefficients at the price 
of extending the so defined loop algebra to all such loops with any winding number 
(with respect to an orientation on the three torus) by making use of the Mandelstam 
identities 

T^[A]Tp[A] = ^[T,o^ + T„o/3-i] (5.25) 

(0) where o is the product in the loop group (with respect to a given base point). 
This Mandelstam identity becomes here just the trigonometrical identity 

cos(q;) cos(/3) = -[cos(a + (3) + cos(a — 13)] (5.26) 

but we refrain from introducing the winding number because we will have to deal with 

higher order polynomials of loop variables anyway. 

Note incidentally that (5.21) implies that the the conjugate variable to T^^,^-. is given 

by 2arth(T,(^)) ! 

Remark : 

It is clear that the Hamilton- Jacobi functional S := cT^r^z) solves all constraints for the 

■^This result implies, for example, that any topology on A/Q which is based on the use of the loop 
variables (traces of the holonomy) displays A/Q as a non-Hausdorff space (|14 ) ! 
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degenerate sector for any constant c. Thus it turns out that the loop variables serve 
as 'good' coordinates only for (part of) the degenerate sector as for the full theory (||T6|). 



We now have to extract 7(2;). One way would be to make use of the so-called area- 
derivative ([0), whose relevant components for the already defined loops become (the 
area-derivative in the (ab)-plane is denoted 6ab) at the point with coordinate z (x,y are 
irrelevant) 

1 1 /•! 

SxyTa{z) = --Ax{z)Ay{z)sm{- dz-f{z)) (5.27) 

SyzT^ix) = -lAy{z)^{z)sm{^Ax{z)) (5.28) 

S..T^iy) = -^iiz)A^iz)sm{^Ayiz)) (5.29) 

but since the area-derivative is a certain limit of a difference between one of the above 
loop variables and an extended one (the area derivative in the (ab)-plane with respect 
to a given loop at one of its points is obtained by attaching an infinitesimal loop lying 
in the (ab)-plane to the given loop at the point of interest, computing the trace of the 
holonomy for this modified loop and the original one, taking the difference and dividing 
by the area element) we actually leave the above algebra and look at an infinite number 
of new loops. We would like to have only as many loop variables as reduced canonical 
variables. Accordingly, we are looking for a better suited object. 
After some trial, one of the most simple objects seems to be an eyeglass loop (0) 
obtained by drawing two loops with equal orientation at zi and Z2 in x or y direction 
connected by a line in z-direction back and forth between zi and Z2. 
We will denote this loop by a{xz) or <7{yz). The computation of the associated T° 
reveals 

11 p2 1 1 

Ta{xz){zi,Z2) = cos{-Ax{zi))cos{-Ax{z2))-cos{ dz-f{z)) sm{-Ax{zi)) sm{-Ax{z2)) 
Z Z J ZI z z 

(5.30) 
and thus allows to extract 7(2;) by taking the derivative with respect to, say, Z2 = z. 
Fix a value Zq on the z-circle and write 

Ta{xz){z) ■■= T^(xz){zo,z) . (5.31) 

We observe that 

T.ixz)izo + 1) = [T^(x.)(^o)]' - (2[T^(.)]' - 1)(1 - [r.(.)(^o)]') (5.32) 

whence 



Ta{z) 



1 — Ta(xz)izo) 

^2(i-[r.(.)(^o)]^) 



(5.33) 



so that the information contained in T„(^z) is also contained in T„(^xz){z) (if one uses 
also To-(^)(2;)) (note that even its sign can be determined from (5.23) as the T^(^a){z) are 
all positive). Therefore we will consider T„{z) as redundant in the following although 
we will use it as an abbreviation for (5.33). Also the analogon T„{yz) oi T„(^xz) is not 
needed. 
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We now compute the extended loop algebra : the only non-vanishing brackets of our 
new element with the rest are (up to a sign) 

{T.(.,)(z),r^%)(z')} = -X[.o,.](^')sin( r7(t)rft)sin(i f' dHt)) 

xy/(l - [T.(.)]2)[(l - [T.(.)(zo)]2)(l - [T.(.)(^)]2) - (T.(..)(z) - T.(.)(^o)T.(.)(^))2] 

(5.34) 

{t.(.,)(z),t;(,)(/)} = i(5(z',zo)[-T;(,)(z) + t;(,)(zo)t.(..)(^)] 

+5(;^', Zi)[-T;(,)(zo) + T;(,)(^)T.(.,)(^)]) (5.35) 

Unfortunately, the right hand side becomes quite complicated and it would seem to 
be better to work with an overcomplete set and apply the algebraic quantization pro- 
gramme (|]l3])- III our case that would boil down to the problem to find loop- variables 



that depend only on the configuration variables, contain sm{^A^{z)),sm{jQ dz^{z)/2) 

and yield brackets that have analytical right hand sides for the corresponding Poisson- 

brackets, but that turns out not to be straightforward without regarding the whole set 

of loops. 

We now come to set up a bijection (up to signs) between the Loop-variables defined so 

far and the canonical variables determined in section 3. 

We have already defined the T's in terms of A^, 11^; 7, 11'^. Here is the local 

inversion : 

A^{z) = 2arcos(T^(^)(2;)) (5.36) 

T'^ (z) 
W'iz) = , "^^^^ (5.37) 

sJl-lT.Miz)? 

j{z) = [arcos( , ^-w(^o)^-w(") - ^-(-)(") )]^ (5.33) 

y'(l-[T.(.)(z)]2)(l-[T.(.)(zo)]^) 

n^W = 7f f ^ , , (5.39) 

^ 2(l-[T,(,){^o)F) 

where the prime denotes derivation with respect to the z-argument only. The z- 
derivation considered as a difference of Loop- variables does not lead out of the class of 
variables that we have defined and so we are fine. 

(5.36)-(5.39) defines a (local) point transformation which is always a local symplecto- 
morphism. Accordingly, when plugging the above formulae into the symplectic poten- 
tial, we are actually able to determine the momenta conjugate to the three variables 

rp rp I 

-'- cr{^) ; -^ cr(xz) ■ 

After some tedious algebraic manipulations, an integration by parts, using that Jgi '-y{z)dz 
written in terms of T°'s is a spatial constant and making use of the fact that quantities 
frozen at zq are also spatial constants we derive the momenta conjugate to 
Ta(x){z),T(j(^y){z),Ta(^xz){z), the basic loop variables that we are using, respectively as 

{T^iz)i^)]'[Ta{x)izo) - T„ix){z)T„ixz)iz)] 



^(1 - [T.(.)]2)[(l - K(.)(zo)]2)(l - [T.ix){zW) - {T.ixz){z) - T^ix){zo)T.ix){z)r 
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X 



1-(T.(.)(.))2 1-[T.(.)(^)]2' ^ • ^ 



T!:u^ 



K(^)(^)]' 



(5.41) 



^(1 - [T.(.)]2)[(l - [T.(.)(2;o)]^)(l - [T.(.)(^)]2) - (T.(..)(2;) - T.(.)(^o)T.(.)(2;))2] 

(5.42) 

and the momentum conjugate to T^(x){zo) is the integral over S^ of expression (5.40) 
only that the roles of T„(^x){zo) and T„(^x-){z) are interchanged. The reason for the ap- 
pearance of this global degree of freedom is that we fixed the point zq to be an observable 
by hand. 

Next we look at the constraints expressed in terms of loop- variables. From their defi- 
nition it is obvious that To-(^) is a scalar, and since (5.38) shows that the density 7 is 
a spatial derivative of a function of our three loop variables, we find that also T„(^x,z) 
must be a scalar (alternatively, this follows from its definition). T„(^x){zo) is a constant 
and its conjugate, as an integral over a density, is also diffeomorphism invariant. 
Thus, necessarily the vector constraint must be 

{T^i^^yS^ + {T^^x,)yS (5.43) 

where we have denoted the corresponding conjugate momenta in (5.40), (5.41) by S'^, S. 
Of course, it is also obvious from the definitions that T^^t . and T^^-j are densities and 
a scalar respectively. 

It turns out to be more convenient to make direct use of the inversion formulas (5.36)- 
(5.39) to obtain the explicit expressions : 

r-pl rpx rpl rpx 

Y _ 2\ o'(^) <^{^) I ^(y) o'{j/) 



'l-{T^(x)y l-{T^(y)] 



2 



^<T{a;)^CT{a;) T^{xz) , {^a{z)) 



+ [arcos(^= ^^"^ " "^-^ ^)V ^ ^^^>' 1 (5.44) 

7(1-[T.(.)F)(1-[T°(,)P)\ yi^^ ^ 



C = 4[ 



arcos(T,(^) )T;(^) arcos(T,(j,) )T^(j^) 



2 



^arc os(T^fa))Tj^(^) arcos(T^(^))T;(^) ^ 



X [arcos(^^^fe±I^i^^^±L)]' ^^^=] (5.45) 

':i-[T.^x)V){i-[T^^x)?) , y^^^ 



^ 1-[T" J2 

L cr(x)J 

where all expressions are to be evaluated at z except for those with a superscript 
which are to be evaluated at zq. 

The result looks much more complicated than when using the variables of the previous 
section. From the point of view of the loop representation one would have to take 
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equations (5.44) and (5.45) and apply it (after a choice of ordering) to a loop state 
(0) according to the Rovelli Smolin loop representation of T", T^. In order to do that 
one would have to expand the trigonometrical formulas involved in both constraints 
and one would obtain an expression involving an infinite sum of loop states. This is a 
very discouraging result, the more as it seems that the expressions (5.44) and (5.45) 
are non-analytical in nature so that it is unclear how to define the action of both con- 
straints on loop states in the loop representation. 

It is clear how to obtain the observables of the theory : simply substitute the loop 
variables via (5.36)-(5.39) into the expressions for the observables given in section 3. 
These formulas are too long and it is not worthwhile as to display them here. 



6 Conclusions 

Let us briefly summarize what we have learnt by studying this model system : 
First of all we worked out a number of ideas and methods to reduce a spatially dif- 
feomorphism invariant field theory which might be useful to keep in mind for future 
work. In particular, our methods do not suffer from the problem of the existence of 
a bijection between a frame on the hypersurface S and a collection of scalar field to 
solve the diffeomorphism constraint as it occurs in the second reference of 0|). 

The symplectic reduction of the diffeomorphism constraint has expectedly lead to 
nonlocal 'already smeared' Dirac observables in a natural way. It was not necessary to 
make use of a complete system of mode functions on the hypersurface (exploiting the 
fact that the Hilbert space of square integrable functions on S is separable) to write 
the true degrees of freedom in such a countable fashion. 

The solution of the scalar constraint quadratic in the momenta reveals that the 
Dirac observables will in general mix momentum and configuration variables relative 
to the polarization of the phase space that one started with. 

The formalism of deparametrization applied to field theories shows that there are 
even more serious problems than those that have been discussed for finite-dimensional 
examples in |Q at least on the technical side. 

The fact that the loop- variables are so badly suited for our model is quite discourag- 
ing and scratches at the hope that many authors had when employing them for gravity, 
namely that they could serve as the natural 'polarization' to solve quantum gravity 
non perturbatively. Actually, the quantum solutions that have been found in [0 are 
solutions for the degenerate sector only. In our model, the variables Tf^(z) and T5^^ 
are simple expressions and the former is easily seen to be a Dirac observable on the 
degenerate sector only. However, it is not the only solution of the constraints (in it 
was speculated that the found solutions are complete). The assumption arises that it 
will be a general feature that the solution of the constraints expressed in terms of the 
loop variables will take a simple form only on the degenerate sector. 
On the other hand, one should never overestimate the results obtained from a model 
and it could in fact also happen that our conclusions are due to an artifact of the 
reduction process. 

Especially the fact that we are dealing with a restricted class of loops seems to con- 
tribute to the complicate appearance of (5.44) and (5.45). 
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A Calculation of the spin-connection 

We will derive here a more general result, namely how the spin-connection looks like 

for the Gowdy models ([§]). 

The densitized triad for the Gowdy models takes the form 

/ E^ Ef \ ( EfEf EfE^j \ 

(Ef) = Ef El so that ^"' = EfE^j E^jEj (A.l) 

\ E' J \ {E'Y J 

is the expression for the twice densitized inverse metric E^E\. We conclude that 

[det(g)]' = [<\ei{E1)f = {E' det{Ef)f =: (E'E)'' (A.2) 

and can invert expression (A.l) to obtain the metric tensor 



E' 



( E^E^j -EfE^j 



E 



{qab) = ^ -EfE^j EfEf . (A.3) 



E \2 



V ,^ 

From this we infer the expressions for the triad and its inverse, that is 

et = ^== and e^ = g,,e^ (A.4) 

det(g) 



namely 

/ Ef Ef \ r^ I El -E- 

(en = -j^f^ ^2 El and (ej,) = \— -El Ef | (A.5) 

^^^ [ E^ J V ^ V Of 

which we have to insert into the formula for the spin-connection 

r^ = -^e'^'e^(2ef,,,] + e'^eU^,) . (A.6) 

Note that e^;, = {eD'6bz in the course of the calculation and that e!,e|e*-''' = e!,(e|)'e*-^'' = 
e!,ef = etc. due to the symmetry properties of the triads so that we can simplify 
(A.6) without reference to the explicit form of the triads to 

K = -le^'^ieliiei)' + e^e[{e[y) + <5..e5(e,^)'] . (A.7) 
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More explicitly 






VI = V(e|(e^)' + eJ(eJ)'). (Ai 



2 

Now we have to go into details and have to make use of the expressions (A.l). The 
calculation is rather tedious. The result is given by 

(K) = ^{^-^^ + ^myE-{EyyE^Ey + {E^,YEyEy], 

^^^ + ^myE + {EfyE^Ej - {EtyEyE%o) 

(r;) = ^{-^-^^ + ^[-{E^,yE-{EiyE^Ey + {EyyE^E^], 



-[-{EiyE + {ED'EfEy - [EyyE^E^], o) 



E^ E' 

n = ^{^,^,Ei{Eyy-Ey{E'ty). (a.9) 

These are the corresponding real parts for the Ashtekar connection for the Gowdy 
models. Let us restrict now (A.9) to our model, that is 

E'^cos{(3) -Eysm{(3) 

E^sin(/5) Eycos{l3) | . (A. 10) 

Q E' 



— r^ sin(/9) rj^cos(/9) 
r,cos(/3) r^sin(/3) I (A.ll) 

or. 



1 EyE"" 1 E^'E'' 

r^ = ^(±^j^)' r = — ^—{t^-E^y r, = -i3' . (A.12) 

2Ey ^ E^ 2E^ ^ Ey 
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